Evolution of a Large Fermi Surface in the Kondo Lattice 



Junya Otsuki^, Hiroaki Kusunose^, and Yoshio Kuramoto^ 
^Department of Physics, Tohoku University, Sendai 980-8578, Japan 
^Department of Physics, Ehime University, Matsuyama 790-8577, Japan 
(Dated: November 17, 2008) 

Single-particle spectrum of the Kondo lattice model is derived with use of the continuous-time 
quantum Monte Carlo method, combined with the dynamical mean-field theory. Crossover behavior 
is traced quantitatively either to a heavy Fermi-liquid state or to a magnetically ordered state from 
the local-moment state at high temperatures. The momentum distribution in the low-temperature 
limit acquires a discontinuity at the location that involves the local-spin degrees of freedom. Even 
without the charge degrees of freedom for local electrons, the excitation spectra exhibit hybridized 
bands similar to those in the Anderson lattice. Temperature dependence in the zero-energy compo- 
nent of the self-energy is crucial in forming the Fermi-liquid state with the large Fermi surface. 

PACS numbers: 75.20.Hr, 71.27.-|-a, 71.10.-w 



A long-standing problem in condensed matter physics 
is how to resolve the dichotomy between localized and 
itinerant characters of electrons in solids. For many met- 
als with / electrons, and also some with d electrons, 
the strong Coulomb repulsion between localized electrons 
suppresses charge fluctuations at each site, leaving only 
a spin and/or an orbital degrees of freedom. These lo- 
calized degrees of freedom interact with conduction elec- 
trons, which are delocalized over the entire crystal. The 
itinerant and localized electrons affect each other, and 
realize rich physical phenomena such as heavy electrons. 

The simplest model to describe the situation is the 
Kondo lattice model (KLM) given by 

H = ^ekcl^Cka + J^Si ■ (Ti. (1) 

k(T i 

Here Si represents the localized spin of the valence elec- 
tron at the i site, and <Ti = J^aa' '4^'^(jcr"^io-' is the spin 
operator of the itinerant conduction electron. The anti- 
ferromagnetic exchange interaction J > under realistic 
condition is much smaller than the bandwidth, and each 
localized spin acts as a weak scatterer for the conduc- 
tion electrons at high temperatures. As temperature de- 
creases, interactions between the localized spins become 
significant, which are mediated by conduction electrons, 
and a magnetic long-range order may be realized. On the 
other hand, the localized spins also tend to be quenched 
by conduction electrons, which is called the Kondo ef- 
fect. If the system remains paramagnetic, coherent quasi- 
particles may emerge by the collective Kondo effect. Al- 
though the overall picture [l[ mentioned above is widely 
accepted, the rich crossover phenomena due to emergent 
quasi-particles remain highly nontrivial ■ 

The KLM has no apparent counterpart of non- 
interacting system, since the perturbation series with re- 
spect to J is essentially singular at J = 0. Accordingly, 
we have no clear starting point to study the nature of the 
paramagnetic ground state in the KLM. On the other 
hand, an ordinary Fermi-liquid argument is applicable 
to the Anderson lattice, since the adiabatically contin- 
ued non-interacting ground state is evident [|. A clue 



for characterizing the KLM is the Luttinger's theorem 
which relates the volume of the Fermi surface with the 
number of electrons j] . A version of proof [5] of this the- 
orem states that the volume of the Fermi-surface is the 
same as that in the Anderson lattice (the so-called "large 
Fermi surface") with unit occupation of each local state. 
The proof is valid provided that the ground state of the 
KLM is a Fermi liquid. However, it is not trivial whether 
this condition is satisfied. In fact, another possibility of 
a "small" Fermi surface has also been proposed d, 0, Q 
with inclusion of explicit coupling between local spins. 

Because of the difficulty in including higher-order ef- 
fects of J, the KLM has been studied either by vari- 
ants of mean- field theory 0,[l[ll, or in one dimensional 
systems [ll|, d, E^l- the mean- field theory, the col- 
lective Kondo effect appears as a phase transition, and 
it is difhcult to reproduce the Kondo crossover at finite 
temperatures. As a particular feature in one dimension, 
on the other hand, the low-energy excitations behave as 
in a Tomonaga-Luttinger liquid without discontinuity in 
the momentum distribution. Furthermore, a magnetic 
long-range order, which competes with the paramagnetic 
state, is strongly suppressed. Thus it is highly desirable 
to obtain reliable information in higher dimensions, es- 
pecially the temperature-dependent evolution of the col- 
lective Kondo effect. 

In this Letter, we report results which become exact 
in the limit of infinite dimensions. First, evolution of 
the Landau quasi-particles is traced quantitatively from 
the local-moment regime at high temperatures. Secondly, 
magnetic phase diagram of the KLM is derived. The pe- 
riodic lattice effect is dealt with by the dynamical mean- 
field theory (DMFT) JA]. To solve the effective impurity 
problem in the DMFT, we use the continuous-time quan- 
tum Monte Carlo (CT^QMC) algorithm adapted 
to the Kondo model [17]. Since no approximation such as 
discretization is involved in the CT-QMC algorithm, and 
since the simulation does not encounter the minus sign 
problem, the effective impurity problem can be solved 
with desired accuracy. 

Model and method. — Let us start with the 
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FIG. 1: The phase diagram of the KLM at ric = 0.9. The in- 
tensity map represents ^ defined by Eq. (O, and the coherent 
energy scale T* is defined by ^ — 0.5. 
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conduction-electron Green function in the KLM: 



(2) 



where e„ = (2n + 1)'!tT is the ferniionic Matsubara fre- 
quency. The self-energy Ec(ie„) takes account of non- 
perturbative contributions from the localized spins, but 
the momentum dependence has been neglected according 
to the DMFT. The KLM is then mapped to the impurity 
Kondo model in the effective medium, which is charac- 
terized by the cavity Green function, 5°(ie„). The im- 
purity self-energy Sc(ie„) should be identical with that 
in Eq. Q , and is given in terms of the impurity t-matrix 
t{ien) as Y.ci'ien)^^ = t{ien)~^ +t/c Here ^c(ie„) is 

determined from the local Green function Gc{itn) of the 



lattice system via tj"(ie„) = Gdie-n 



Sc(«e„). 



We adopt a nearest-neighbor tight-binding model in 
the infinite-dimensional hyper-cubic lattice ^18] , which is 
characterized by the Gaussian density of states: Pc(^) = 
D-^y/2/TTexp{~2Lu^/D^). We take D = 1 as the unit of 
energy, and fix the conduction-electron density per site 
as He = 0.9 throughout this paper. The latter choice fa- 
vors antiferromagnetism by the nearly nesting condition. 
We use 10^ QMC samples at most with 10-100 intervals. 
The DMFT self-consistent equations converge typically 
within 10 iterations. 

Phase diagram. — We first show the J-T phase 
diagram in Fig.[T] The antiferromagnetic transition tem- 
perature Tap is determined by the divergence of the stag- 
gered magnetic susceptibility, which can be computed in 
terms of the Bethe-Salpeter equation with numerically 
derived local vertices [3, [l^. Here we have neglected 
possible incommensurate order for simplicity, since our 
main interest in this paper lies in the paramagnetic phase. 
The overall structure of the phase diagram is understood 
by competition between the Kondo effect and the RKKY 
interaction Q. As J increases (T decreases) in the para- 
magnetic phase, the localized moments gradually disap- 
pear, and the Fermi-liquid behavior dominates below the 
coherent energy scale T*. The intensity map in Fig. [1] 
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FIG. 2: The single-particle excitation spectrum A{k, ui) for 
J = 0.3 and Uc = 0.9 at (a) T = 0.25 and (b) T = 0.0025. The 
slanted line represents the non-interacting spectrum ui = K — jj. 
which is realized with J = 0. 



represents the degree of itinerancy in terms of the pa- 
rameter ^, whose definition is given later in Eq. ([6]). We 
will show that the Fermi liquid is indeed realized for 
J > Jc ~ 0.2 at low temperatures. The crossover re- 
gion in temperature becomes narrower near the quantum 
critical point, Jc- 

Single-particle spectra. — In the DMFT, k- 
dependence of the Green function enters only through e^. 
Therefore Gc{k,ien) is regarded as a function oi k — Ck 
and is written as 



(3) 



The single-particle excitation spectrum j4(k, oj) = 
— ImG'c(K, ijj+iQ)/!: is obtained by means of the Fade ap- 
proximation [20| . The validity of the Fade approximation 
has been confirmed for the impurity Kondo-type mod- 
els • Our QMC data has the accuracy high enough to 
obtain reliable real-frequency spectra. 

Figure [2] shows the single-particle excitation spectra 
for a fixed value of J = 0.3. We define two energies, 
kl and ksi corresponding to the large and small Fermi 
surfaces, respectively. The upper panel shows A[k, uj) for 
T = 0.25. The spectrum exhibits a behavior of almost 
non-interacting electrons at high energies. However, the 
coupling with spin fluctuations give rise to a pseudo-gap 
near the small Fermi surface ks- 
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The lower panel shows the spectrum at T = 0.0025, 
which is much lower than the impurity Kondo tempera- 
ture defined by Tr = ^/?e-l/9 - 0.1 with g = 2Jpc(0). 
The pseudo-gap looks like a true gap, and the coherent 
quasi-particles develop near the gap edges. This struc- 
ture has a close resemblance to the hybridized band in the 
Anderson lattice. In particular, the width of the pseudo- 
gap is the order of the Kondo temperature Tk ~ 0.1. The 
quasi-particle intensity becomes sharper near the Fermi 
energy, and the quasi-particle band crosses the Fermi en- 
ergy exactly at kl. Namely, the localized spins do con- 
tribute to the Fermi volume. 

Evolution of the Fermi liquid. — Let us charac- 
terize the Green function near the Fermi surface. We 
expand the self-energy in powers of ie„ to obtain 



Gc(K,ie„) 



ien - z{k - m + Sc(0)) 
where the renormalization factor z is given by 



z= 1 



9ImI]c(ie„) 



der, 



(4) 



(5) 
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From the QMC data, we have confirmed numerically 
the Fermi-liquid conditions, 9ReI]c(?e„)/c)e„ and 
Iml]c(ie„) ^ in the limit e„ — > 0, which means T — > 0. 
As a result, we obtain z ~ 0.084 for J = 0.3. We extrap- 
olate ReSc(O) from the QMC data of ReEc(ien) for small 
|n|, assuming a quadratic function centered at e„ = 0. 

We now demonstrate the evolution of the Fermi-liquid 
state with the large Fermi surface as a function of tem- 
perature. The Fermi surface appears at such momentum 
that satisfies Cfc = /x — Ec(0). Figure [3] shows temperature 
dependences oi fi — ReSc(O) and fi for J = 0.2, 0.3 and 
0.4. There are apparently two regimes in temperature. 
At high temperature, the self-energy correction is not im- 
portant, and thermal excitations occur around /i ^ ks- 
At low temperature, on the other hand, /.t— Rel]c(0) devi- 
ates significantly from fi, and gives rise to the large Fermi 
surface at kl. From this observation, /i — Rel]c(0) turns 
out to be a good measure to characterize the crossover 
between the local-moment regime and the Fermi-liquid 
state. We introduce a dimensionless quantity 



^ = (/i - Rel]c(0) - kl)/(ks - 



(6) 



which ranges from to 1 for T <C The limit ^ = 
represents the Fermi-liquid state with the large Fermi 
surface, and the other limit ^ = 1 represents almost non- 
interacting conduction electrons. We define the crossover 
temperature T* by the condition ^ — 0.5 as shown in 
Fig. [H At J = 0.3, for example, we obtain T* ~ 0.035, 
which is smaller than Tk ~ 0.1. 

For the one-dimensional KLM, the existence of two 
different energy scales has already been discussed in the 
literature [2l| . The uniform charge susceptibility Xc{T), 
for example, takes maximum around Tk and decreases 
with decreasing temperature due to the development of 



w 

I 




FIG. 3: Temperature dependences of the chemical potential 
fj, and its self-energy correction, /x — ReEc(O), for ric = 0.9. 



the pseudo-gap. With a further decrease of temperature, 
XciT) turns to increase reflecting contribution of coher- 
ent excitations, and it saturates to a finite value. For 
the infinite-dimensional KLM, our results of the uniform 
charge susceptibility, which will be reported elsewhere, 
shows similar behavior with that in one dimension. The 
distinct energy scales without phase transitions are im- 
portant to analyze experimental data properly. 

We note the difference between the present descrip- 
tion in Eq. ^ and the ordinary Fermi-liquid state in the 
Anderson lattice In the Anderson lattice, the local- 
ized electrons which are practically immobile at high tem- 
peratures eventually become itinerant. Thus, the major 
component of the quasi-particle is the localized electron. 
On the contrary, the local spins in the KLM remain ab- 
solutely immobile although the Fermi surface becomes 
large. At high temperatures, the conduction electrons 
interact incoherently with local spin fluctuations. At 
low temperatures and energies, the conduction electrons 
manage to form coherent quasi-particles which accom- 
pany the polarization cloud of spin fluctuations. These 
quasi-particles exhibit the heavy-fermion behavior. In 
particular, the origin of the large linear specific heat is 
the entropy of the local spins as in the Anderson lattice. 
The spin entropy changes in the temperature scale of Tk. 

Momentum distribution. — We now discuss the 
momentum distribution of conduction electrons: 

(cLcfe.) = T^Gc(«;,ze„)e^^"°+ = n,{K). (7) 

n 

Figure |4] shows nc{K) in the Fermi- liquid state, T <^T* . 
The gradient of nc{n) at k = kl becomes steep with 
decreasing temperature, indicating the existence of dis- 
continuity at T = 0. The magnitude of the discontinuity 
is consistent with z ~ 0.084 estimated from Ec(ien). On 
the other hand, nc(K) at ks is almost independent of T 
in the Fermi-liquid state. For characterization of the mo- 
mentum distribution at finite temperature, we introduce 
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FIG. 4: Momentum distribution nc(K) for J = 0.3 and = 
0.9. The vicinity of the large Fermi surface is enlarged in the 
left inset. The right inset shows the temperature dependences 
of the "width" of nc(K) at k = res and kl. 



the energy scales (7 — L, S): 



t; = - 



Ok 



(8) 



at kl and Kg . In the Fermi- liquid state with the approxi- 
mate Green function, Eq. (jH), Tl should behave as T/z^. 
The inset in Fig. |3] shows the temperature dependences 
of T-y. The linear-T dependence in Tl ~ T/z^ at tem- 
peratures lower than T* (~ 0.035 for J — 0.3) verifies 



the existence of the discontinuity at T = 0, and ensures 
the validity of the quasi-particle expansion in Eq. ([4|. 
At temperatures higher than Tk, on the other hand, kl 
loses its importance. Instead, thermal excitations are 
populated around /i ~ ks which develops according to 
Ts~T. 

We have thus established the Fermi-liquid state of the 
KLM in infinite dimensions. On the basis of the Lut- 
tinger's theorem, we may further develop a phenomeno- 
logical Fermi-liquid theory in the vicinity of the large 
Fermi surface. Note that the Fermi-liquid state realized 
in the KLM has no explicit counterpart of non-interacting 
system. This situation is identical to the local Fermi- 
liquid theory by Nozieres in the single-impurity Kondo 
model [12] , where the phase shift of conduction electrons 
is a key quantity. Following similar strategy, an extension 
should be possible starting from the large Fermi surface 
in the three-dimensional KLM. Investigation of possible 
superconductivity is also an challenging issue in extend- 
ing the present infinite-dimensional theory for the KLM. 

In summary, we have demonstrated the temperature- 
dependent evolution of the Fermi-liquid state with the 
large Fermi surface. Our numerical results in the low- 
temperature limit confirm the premise for Luttinger's 
theorem that relies on the Fermi-liquid character of low- 
energy excitations Q. 
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